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Abstract 

We study general mathematical framework for variation of potential energy with respect 
to domain deformation. It enables rigorous derivation of the integral formulas for the energy 
release rate in crack problems. Applying a technique of the shape sensitivity analysis, we 
formulate the shape derivative of potential energy as a variational problem with a parameter. 
Key tools of our abstract theory are a new parameter variational principle and the classical 
implicit function theorem in Banach spaces. 
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1 Introduction 

Many variational problems related to the domain deformation have been investigated in the 
theory of shape derivatives or the shape sensitivity analysis ([9], [E]). According to the spread 
of the importance of shape derivatives in various scientific fields, more development of their 
mathematical foundation has been required. The purpose of this paper is to establish some 
abstract parameter variational formulas and their application to the shape derivative of potential 
energies. An important example is the energy release rate in crack problems, which is known as 
one of the most fundamental quantities in the theory of fracture mechanics. 

Scientific investigation to understand crack evolution process in elastic body was originated 
by Griffith [6] and has been studied from various viewpoints in engineering, physics and mathe- 
matics since then. Griffith's idea in the fracture mechanics is even now the fundamental theory 
in modeling and analysis of the crack behavior. We here make reference to only very few ex- 
tended studies from mathematical point of view, Cherepanov [3], Rice |14j . Ohtsuka [10], 
[12j , Ohtsuka-Khludnev [13] , and Francfort-Marigo [5] . For more complete list of crack problems 
and fracture mechanics, please see the references in the above papers. 

In the Griffith's theory and its various extended theories such as [5], the concept of the energy 
release rate G plays an important role. According to such theories, we treat crack evolutions in 
brittle materials with linear elasticity under a quasi-static situation, in which applied boundary 
loading is supposed to change slowly and any inertial effect can be ignored. The elastic energy at 
a fixed moment is supposed to be given by minimization of an elastostatic energy. According to 
the Griffith's theory, the surface energy required in the crack evolution is supplied by relaxation 
of the potential energy along crack growth. 

Roughly speaking, the energy release rate G is defined as follows. Please see the above 
references for more precise definition. Let f2* be a bounded domain in R n (n > 2), which 
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corresponds to the uncracked material under consideration. We assume that a crack E exist in 
where E is the closure of an n — 1 dimensional hypersurface. The cracked elastic body is 

represented by f2* \ E. We consider a virtual crack extension E(i) with parameter t £ [0,T), 

where E = S(0) C E(ti) C E(i 2 ) (0 < v t x < v t 2 < T). 

Under the quasi-static assumption, the elastic potential energy E(t) in £l(t) := \ E(t) is 

given by 



E(t) := min / t> (x), Vv(x))dx, (1.1) 

11 ifi,\S(i) 

where by W(x, f (x), Vti(i)) the potential energy density including a body force is denoted, and 
min„ is taken over all possible displacement fields in $7* \ E(i) with a given boundary condition. 
For the admissible displacement fields, a given displacement field is imposed only on the part 

C On the other part d£l(t) \ including both sides of E(i), the normal stress free 

condition is imposed for the minimizer on d£l(t) implicitly. 

The energy release rate G at t = along the virtual crack extension {E(t)}o<t<T is given by 

G := lim ?°»- g( f ) . (1.2) 

i^+o |E(t)\E| v ; 

Since E(t) < E(0), G > follows if the limit exists. The Griffith's criterion for the brittle crack 
extension is given by G > G c , where G c is an energy required to create new crack per unit length 
and it is a constant depending on the material property and the position. 

Cherepanov [3] and Rice [TT] studied so-called J-integral for straight crack in two dimensional 
linear elasticity, which is a path-independent integral expressions of the energy release rate. Since 
these works, theoretical and practical studies of crack evolutions have been much developed by 
means of such useful mathematical expression of G in two dimensional case. 

As an alternative approach to such energy based arguments, Irwin [8] proposed the notions 
of fracture toughness and stress intensity factors and he developed arguments based on the 
singularity of stress fields. 

While most of these mathematically rigorous results have been restricted to two dimen- 
sional linear elasticity (and often only for straight cracks), Ohtsuka [JU], [TT], [12] and Ohtsuka- 
Khludnev [13] developed a mathematical formulation of the energy release rate for general 
curved cracks in multi-dimensional linear or semi- linear elliptic systems. They proved existence 
of the energy release rate, and obtained its expression by a domain integral and by a generalized 
J-integral. 

Based on the idea in [TT], we shall give a new mathematical framework for shape derivative of 
potential energy including the energy release rate. Adopting domain perturbation (p of Lipschitz 
class, we treat the shape derivative as an abstract parameter variation problem in Banach spaces, 
where (p is considered as a parameter belonging to a Lipschitz class. Instead of estimating the 
limit (II. 2p directly as in [TT] and |13| . we treat it by means of the Frechet derivative. 

In our approach, the shape derivative of minimum potential energy is derived as a Frechet 
derivative in a Banach space within an abstract parameter variation formulas and it is given as 
a domain integral. The key tools in the abstract parameter variation setting are the implicit 
function theorem and the Lax-Milgram theorem. 

The organization of this paper is as follows. Abstract parameter variation formulas are es- 
tablished based on the implicit function theorem in Banach spaces in Section [2l In Section [3l a 
framework of Lipschitz deformation of domains, which includes crack extensions, is introduced. 
Minimization problems with a general potential energy in deformed domains are studied in Sec- 
tion H] as an application of the abstract parameter variation formulas in Section [2j Quadratic 
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energy functionals corresponding to second order linear elliptic equations are treated in Sec- 
tion [5j Under a weak regularity assumption, we show the existence of the shape derivative of 
the minimum potential energy and derive its domain integral expression and a boundary integral 
formula ( J-integral) . The results obtained there include the results in [11] and [13] under a weaker 
assumption for regularity of domain perturbation. In [13j . they assumed that the domain pertur- 
bation t — » ip(t) belongs to C 2 ([0, T], W 2,co (W l ) n ) and derived the domain integral expression, 
whereas, in Theorem 15.51 we prove it under a weaker assumption ip E C 1 ([0, T], W 1 ' 00 (M n ) n ). 
For simplicity, description will be made on scalar equation in Sections H] and [5j But our results 
are easily extended to elliptic systems such as linear elasticity problems as shown in [llj. 

2 Parameter variation formulas 

We consider a variational problem with a parameter in an abstract setting. For a real valued 
functional J defined on a metric space S, uq E S is called a global minimizer of J in S, if 
J(uq) < J(u) for all u E S. If there exists an open set O C S and uq is a global minimizer of J 
in 0, no E S is called a Zoca/ minimizer of J. 

Let X and M be real Banach spaces. For open subsets Uq C X and Oq C M, we consider 
J E x Co 5 lR) and n E C 1 (Oo ) ^o)- We assume that n(n) is a local minimizer of J(-,u) in 

for each u E 0o> and define J*(n) := ./(«(«),«) for n E Oq. Then we have J* E C 1 (Oo) and 



J*(n) = £yj(n(n),n)] = d x J(u(n),n)[u'(n)] + d M J(u{y),n) = d M J(u(n),n), (2.1) 



where «/£ denotes the Frechet derivative of J* and denotes the Frechet differential operator 
with respect to u E M. The symbols dx and 8m denote the partial Frechet derivative operators 
for J(n, n) with respect to n E X and n E M, respectively. The last equality of (12. ip follows 
from dxJ(u(fi),fj,) = E X' , where X' denotes the dual space of X. The formula 



is a simple but essential equation in this paper. 

The following fundamental theorem states that the formula (|2.2p is derived under a weaker 
assumption for regularity. 

Theorem 2.1. Let X and M be real Banach spaces. ForlA® C X and an open subset Oq C M, 
we consider a real valued functional J : Uq x Oq — > R and a map u : Oq ^ Uq. We define 
J*(n) := J(n(n), u) for a E Oq. We suppose the following conditions. 

1. J E C°(U x O ), J(w, •) E C^Co) /or w E W , and d M J E C°{U x O ,M'). 

2. n E C°(Oo ; X) w (a0 i s a global minimizer of «/(•,«) in Uq for each u E Oq. 
Then J* E C^Oq) and {UP Zioids. 

Proof. We fix no G Co and we define no := n(no) and 



4(n) = d M J(u(n), fJ .) (nE0 o ) 



(2.2) 



r(n) := J*(n) - J*(n ) - <9 A / J(n , u )[n - no] (« E O ). 



Since n(n) is a global minimizer and n E C°(Oq, X), if a is close to no, we have 



r(n) < J(n , n) - J (n , n ) - <9a/ J (n , n ) [u - u ] 



o(||« - «o||m), 
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r(n) > J{u{n),n) - J(u(n),fj, ) - d M J(u ,fi )[n - fm] 
1 

{dMJ(u(fj.),fio + s(jjl - fio)) - dMJ(uo,no)}\p, - no]ds = o(\\/jl - ho\\m)- 

It follows that r(n) = o(||/x — /^oIIm) as \i — > fio, and we obtain the formula (|2.2[) and J* G 
C°(O ,M'). □ 

Corollary 2.2. Under the condition of Theorem \2.1\ we assume that Uq is open. If 8m J G 
C k {U x O , M') and u G C fc (O ,^), tten J* G C fc+1 (O )- 

Proof. This immediately follows from the formula (|2.2p . □ 
We apply the implicit function theorem in Banach spaces below. The proof is found in [2] 
and [7] etc. For two Banach spaces X and Y, B(X,Y) denotes the Banach space which consists 
of all bounded linear operators from X to Y. 

Theorem 2.3 (Implicit function theorem). Let X, Y , Z be real Banach spaces and U, V be 
open sets in X and Y, respectively. We suppose that F : U x V — > Z and (xo,yo) G U x V 
satisfy the conditions; 

1. F(x ,y ) = 0. 

2. F G C°(U x V,Z). 

3. F(x,-) G C l (V,Z) for x G 17 and By F is continuous at (x,y) = (xo,yo). 

4. (^FfxcJA))) -1 €S(Z,y). 

T/ien i/iere exist a convex open neighborhood of (xo,yo), Uq X Vo C U x V and f G C°(Uo, Vq), 
such that, for (x, y) G Uq x Vo, y) = if and only ify = f(x) . Moreover, if F G C k (U x V, Z) 
(k G N), then f G C fc (C/ ,^o)- 

From Theorem 12.11 and the implicit function theorem, we get the following theorems. 

Theorem 2.4. Let X and M be real Banach spaces andU and O be open subsets of X and M, 
respectively. We consider a real valued functional J :IA x O ^ R and fix fiQ G O. We assume 

1. J(-,n) G C 2 {U) for neO and d x J G C°(U x 0,X'). 

2. uq €U satisfies 8xJ(uq, fio) = 0. 

3. d\,J is continuous at (u>,/u) = (uo,A*o)- 

4. There exists a > such that d\ J{uq, Ho)[w, w] > a\\w\\x for w G X. 

Then there exist a convex open neighborhood of {uq,Hq) ; Uq x Oq C U x O and u G C°(Oq,Uq), 
such that, for [i G Oq, the following three conditions are equivalent. 

a. w G Uq is a local minimizer of </(•,//) 

b. w G Uq satisfies dxJ(w,fi) = 0. 

c. w = u(n). 

In this case, u(fi) is a global minimizer of J(-,fi) on Uq. 
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Proof. We define a map F := dxJ from U x O to X' and apply Theorem 12 .31 at (w, p) = (uo, po). 
From assumption 4 and the Lax-Milgram theorem, dxF(uo, /1q) = d\ J(uq, /jLq) becomes a linear 
topological isomorphism from X to X' . Then, from the implicit function theorem there exist 
a convex open neighborhood of (uq,/j,q), Uq x Oq C U x O and u G C (Oo,Uq), such that, for 
H G Oq , w G Uq satisfies dxJ(w, fi) = if and only if w = u{p). 

From the continuity of d\ J at (uo,po), without loss of generality, (after replacing Uq and 
Co with smaller ones if needed) we can assume that 

d x J(v,n)[w,w]>^\\w\\ x (^el, VrfeWoxOo). (2.3) 

For fx G Co, if if G Z^o is a local minimizer of J(-,p) in Z^o, the dxJ(w,p) = follows. 
Conversely, if w G Z^o satisfies dxJ(w,p) = 0, w is a local minimizer in Uq from the condition 
(|2.3p . It also follows from (|2.3p that u(p) is a global minimizer of J(-,p) in ZYq- □ 



Theorem 2.5. Under the condition of Theorem \2.4\ we additionally assume that dxJ G C k (Ux 
O, X') for some k G N. T/ien u G C k (0 ,Uo) holds. 

Proof. The assertion follows from the implicit function theorem. □ 
Under the condition of Theorem 12.41 we define 



J*(//) := J(u(n),fi) (fj,eO ). 

As a consequence of Theorem 12.51 a sufficient condition for J* G C 1 (Co) is J G C (W x 0) and 
<9x</ G C 1 (Zi x G,X'). However, the condition dxJ G C X {U x 0, X') is not necessary due to 
Theorem 12.11 

Theorem 2.6. Under the condition of Theorem \2.4\ we additionally assume that J G C k (U x O) 
for some k G N, then J* G C k (Oo) and it satisfies (|2.2p . 

Proo/. From Theorem EJJ J* G C x (Co) and ([22]) immediately follows. Since u G C fc_1 (O ,X) 
follows from Theorem 12.51 J* G C k (Oo) is obtained from the formula (|2.2p . □ 
Let us consider the case = 1 in Theorem 12.6^ where J* G C 1 (Oo) is derived under the 
conditions J G C l {U x 0) and J(-,/i) G C 2 {U). In this case, u G C o (0 o ,^o) holds from 
Theorem 12.41 but u G" C 1 (Oq,Uq) in general. In order to obtain u G C 1 (Oo>^o)i we need 
to assume <9xJ G C l {U x 0,X') (Theorem I2.5p . We have Holder regularity of u under the 
condition of Theorem 12.61 with k = 1. 



Proposition 2.7. Under the condition of Theorem \2.4\ we additionally assume that J G C 1 (Ux 
O), then we have 

IK/-0 - uq\\ x = o - A*o|ll/ 2 ) as - /j> \\m —> 0. 
Proof. From the proof of Theorem 12.3^ there exists C > such that 

\\u{fi)-u \\x <C\\dxJ{u ,^\\x' (/i G Oq). (2.4) 

Let po > with {v G X; \\v — uq\\x < Po} C Z^o- For h G X with ||/i||x = 1, ^ £ Oo an d 
p G (0, po]) we have 

J(u + p/i, p) = J{uq,h) + p3x J(^o, aOM + P 2 / (! ~~ s)d 2 x J{uQ + sph,fj,)[h, h]ds. 

Jo 
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d x J(u ,p)[h] 
= d x J(u ,p)[h] - d x J(u ,p )[h] 

p~ 1 (J(u + ph,p) - J(u ,fj,)) - p I (1 - s)d x J{uQ + sph, p)[h,h]ds 



p 1 (J(n + ph,p ) - J(u ,p )) - p I (1 - s)d x J(u + sph, p )[h,h]ds 
P^ 1 I {9mJ(u + ph, po + t(p - po)) ~ d M J(u ,Po + t(p - po))}[p - po]dt 



-p / (1 - s){d 2 x J(uQ + sph, p) - d x J(u + sph, p )}[h,h]ds 
Jo 

For r > 0, we define 

S(r) := {(w,X) G X x M; \\w - u \\ x < r, \\X - p \\ M < r 2 } , 

u(r) := sup \\dMJ(w,\) - dMJ(uo,po)\\M' 

+ sup \\d x J(w, A) — d x J(uo,po)\\g 2 r Xt gj. 

(w,X)eS(r) 

1 /2 

We remark that u>(r) — > as r — > +0. Choosing p := \\p — ^oIIm » we orj tain 
\\d x J(u , p)\\ X / < 2u(p)p (peO , \\p- po\\m < pi). 
Hence, from (|2.4p . we have 

\\u(p) - u \\ x < C\\d x J(u , fjt)\\x' < 2Cuj(p)p = o(p). 

□ 

Under the conditions of Theorem 12.41 d x J{u{p), p) can be regarded as a linear topological 
isomorphism from X to X' from the Lax-Milgram theorem. Therefore, we can define A(p) € 
B(X',X) which satisfies 

d x J{u(p),p)[k{p)h,w] = h[w] ( v ™ eX, v /i G X'). 

The Frechet derivative of the local minimizer u{p) with respect to parameter p is given by the 
next proposition. 

Proposition 2.8. Under the condition of Theorem \2. 51 with k = l, 

u'{p) = -A(p)h (p) (p G O ), (2.5) 

holds, where ho(p) := dMd x J(u(p), p) G B(M,X'). 

Proof. Differentiating d x J(u(p), p) = G X' by p, we have 

d 2 x J{u(p),p)[u'{p)] + dMd x J{u(p),p) = G B(M,X'). 

This is equivalent to (|2.5p from the Lax-Milgram theorem. □ 

Proposition 2.9. Under the condition of Theorem \2.4\ we additionally assume that J G C 2 (Ux 
O) then J* G C 2 {Oq) and it satisfies 

J"(p)[pi,P2] =d 2 M J(u{p),p)[p 1 ,p 2 \- x {A(p)h (p)[p 1 ], h (p)[p 2 ]) X >, 

for p G Oq and p\, P2 G M . 

Proof. Differentiating the formula (j2.2[) by p and substituting (|2.5j) . we obtain the formula. □ 
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3 Lipschitz deformation of domains 

We study a domain deformation with Lipschitz transform 99 : Q — > </?(fi), where Q is a bounded 
domain in M ra (n G N) and <p is a M™-valued Lipschitz function. The identity map on M n is 
denoted by ipo(x) = x (x G 1"). 

For a function u : $7 — > we define 

1 1 - u(y)\ 
\u\up,n ■= sup 1 _ I , 

where | • | denotes the Euclidean norm in R n or M. k . If |«|Lipn < u is called uniformly 
Lipschitz continuous on Q. It is known that, for u G W 1,co (Q), there is u G C°(fl) such that 
u(x) = u{x) a.e. x G SI, in other words, we can regard W 1,oc (Q) C C°(0). If $7 is convex, 
W 1 ' 00 ^) = C ' 1 ^) as a subset of C°(Q). Moreover, if fc = 1, we have 

||Vu|Uoo (n) = |«|Li P ,n (u g w 1 ' 00 ^) n c°(n)). 

In the following argument, we fix a bounded convex domain 0,q C W 1 (n > 2), and we identify 
W 1 ' oo (n ,R n ) with C ' 1 ^,^™). 

Proposition 3.1. Suppose that if G W 1 ' oc (Qq,W 1 ) satisfies \<p — <fo\Lip,ci < 1- Then (p is a 
bi-Lipschitz transform from to (p(Tlo), i.e. ip is bijective from f^o onto an open set and <p and 
(p^ 1 are both uniformly Lipschitz continuous. Moreover, we have 

ess-inf (det Vip T ) > (1 - \cp - (^bpA))" > °- ( 3 - 1 ) 
where V(/? T is the Jacobian matrix defined by 

V<p T (x) := (^(x)) GK nxn (. = (xi,-,x n ) T GO ). 

Proof. Let fi := 92 — 920 and 6* := |/x|Lip,n G (0> !)■ First, we show that 92(^0) is open. We 
arbitrarily fix xo G f2o and define yo := tp(%o)- Let <5 > such that Bg(xo) C J)o ; where 
Bs(xq) := {x G M n ; \x — xq\ < 6}. For y G B( 1 _ } S (yo), we show that y G 92(^0)- It is easily 
checked that := y — /•*(£) is a uniform contraction on Bg(xo). From the contraction mapping 
theorem, there is a fixed point x = T(x) = y — fi(x) in B$(xq), that is y = <p(x). Hence, ip(Qo) is 
an open set. Since \ip(xi) — 92^2) | > \%i — ^2 1 — IM^i) — mO 3 ^)! > (1 — — £2! for xi, X2 G Oo, 
it follows that 99 is injective and satisfies uniform Lipschitz condition on 92(^0). 

From Rademacher's theorem (see [I], |16|). \x is differentiable almost everywhere and the 
derivative coincides with the distributional derivative almost everywhere, i.e., there exists M C 
ft with £ n (M) = such that 

V T /,(x)y = lim & + ^ " (x G tt \ Af), 

h^O h 

where V T /x = (V^ T ) T . It follows that 

|V T /x(x)y| <9\y\ (xen \M), 

and that the moduli of all the eigenvalues of V T /i(x) for x G f^o \ N are bounded by 6. Hence, 
we obtain 

det (V T v?(x)) = det (/ + V T n{x)) > (1 - 6) n (x G n \ AT). 
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_ □ 

We fix an open set Q which satisfies Q C fio an d the deformed domain is denoted 

by f2(</?) under the condition of Proposition 13.14 hereafter. We define a pushforward operator 
ip* which transforms a function d on to a function ip*v := v o tp^ 1 on Q(ip), if <£> satisfies 
Proposition 13.11 We define 



A{fp) := (V^)- 1 G L°°(ft , 



K(<p) := detVv? T G L°°(Qo,] 



These Jacobi matrices and Jacobian appear in the pullback of differentiation and integration on 
Q(<p) to Q. For a function v on Q, we have 



[V(ip*v)]o<p = A(<p)Vv a.einft (« G W 1 ' 1 ^)), 



n(<p) 



v(x)n(ip)(x)dx (v G 



(3.2) 
(3.3) 



These equalities are well known in the case (p G C . However, for </j G C ' 1 , these are not so 
trivial. See, [4] and [16] etc. for details. We omit the proof of the next proposition since it is 
clear from (13.21) and (13.31). 



Proposition 3.2. Under the condition of Proposition HOI for p G [l,oo], (p* is a linear topologi- 
cal isomorphism from L p (0) onto L p (Q((p)), and a linear topological isomorphism from W 1,p {£l) 
onto W lj, (£l((p)). 

The following theorem plays an essential role in the application to the shape derivatives. 
Theorem 3.3. Let $7 be an open subset o/$7o- 

1. k G C 00 (W /1 ' 00 (ri, M n ), -L°°(il)) ; and K'(if )[fA = div/i for fi G W 1, °°(Q,W l ). 

2. AeC 00 (Oo(fi),I 00 (fi,R ,lxn )), where 

Oo(Q) := {9? G W 1,0 °(fi,R n ); ess-inf «(^) > 0}. 

In particular, A'((p )[fi] = -V// T /io/ds /or ^ G W 1,0 °(n, M n ). 
Proof. Since the determinant is a polynomial of degree n, it is clear that k belongs to C°°(W 1,00 (£l, M n ), L°°(J))). 



For fixed G W^°°(n 



we define 



mi<(t) := 5ij + t-f*- G L°°(ft) (i, j = 1, 



where oV,- is the Kronecker's delta. Then we have 



, n, t G 



«'(¥>o)M 



d 
dt 

d_ 

dt 



K(ip + tlx) 



t=0 



d 
dt 



det (mij(t)) 



t=o 



t=o 



■ m r 



t=0 



(m n (t) ■ ■ ■ m nn (t)) = ^2 m n(°) • • • "4(0) • • • m nn{0) = div/i. 



i=l 



Let the (i,j) component of A{ip) be denoted by aij((p) G L°°(Q). Then we have aij{ip) = 
(ip) j «(<£>)) where ctij((p) is the (i, j) cofactor of V(^ T , which is a polynomial of |^ of degree 
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n - 1. Since ess- inf n n((p) > for (p G O (fi), ay G C°°(Oo(0), L°°(fi)) follows. For fixed 
G Ty 1 ' oo (0,M n ), differentiating the identity 

A((po + + tV/~t T ) = / (I: identity matrix of degree n), 

by t G M. at i = 0, we have 

A'{wM+A{ipv)Vf = 0. 

Since A(ipo) = I, we have A'((po)[fj] = — V/x T . □ 
We define an open subset of W 1,00 (Qo, R n ) as 

0(0) := G PF 1 ' 00 ^^); |y> - <A)|LipA, < 1. C ^o} • (3.4) 

Proposition 3.4. VKe assume /j G Vt /1,00 (r2o) wi/i supp(/i) C f2. 

1. // |^|Lip,Q < 1' then ip = (fo + fi is a bi-Lipschitz transform from f2 onto itself. 

2. For t G R mi/* |t/-i|Lip,Qo < 1) u;e define a bi-Lipschitz transform ip(t) = f>o + tfj, from 
Q, to itself. Let I G {0,1} and p G [l,oo]. Suppose that f G W l,p (Q) if p G [l,oo), and 
/ G C z (0) n W 1 ' 00 ^) i/p = oo. Then ip(t)J -> / strono/y in V^' p (ft) as t -> 0. 

Proof. From Proposition [3TT] claim 1 is clear. For claim 2, let us fix to > with |io/^|Lip,Oo 

< 1. 

Then, from Proposition 13.21 there exist C > such that the following inequalities hold for 

|*| < t , 

\\<p(t)*f ~ f\\w''P(Q) = \\<p(t)*(f ~ f y(*))llwi.p(n) < Cll/ - / ° <p(t)\\wi*(ny 
Since [tp ^ / o <p] G C°(0(n), V^< p (ft)), we obtain 

\\f - f O (p(t)\\ W l, p{n) = \\f O(f - f Oip(t)\\ W L p(n) -»0, 

as i -> 0. □ 
4 Potential energy in deformed domains 

In this section, we consider minimization problems of an abstract potential energy in deformed 
domains. Some concrete examples in linear elliptic equations will be given in Section 

Let S7o be a fixed bounded convex open set of W 1 (n G N). We consider an open set £1 
whose closure is contained in fio- For v G H l (Cl) = W 1,2 (Q), we introduce the following energy 
functional: 

E(v,n) := / W(x,v(x),Vv(x))dx, 
Jn 

where 

W(£, 77, C) G M for (£, i),()efioxlx M n , 

is a given energy density function. We assume some suitable regularity conditions and bound- 
edness of its derivatives in the following argument. For simplicity, the partial derivatives of 
W with respect to £, rj and C will be denoted by V{W = (fg,--- , fg) T , W v = f^, and 
V ' (W = (f^r - ' >fp") T ) respectively. Moreover, for v G -ff 1 (0), we often write W(v(x)) = 
W(x,v(x),Vv(x)), v]W(v(x)) = V^W(x,v(x),Vv(x)), etc. 
We consider the following minimization problem. 
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Problem 4.1. Let V be a closed subspace of H l (£l) with Hq(Q) cVcH 1 (Q) ! and let V(g) := 
{v G v — g G V} for g G H 1 ^). For given g G H (fi), /fnd a Zoca/ minimizer u of 

E(-,Q) in V(g), i.e. u G V(g) and there exists p > such that 

E(u,n) < E(w,Q.) Cw G V(g) with \\w -u\\ H i (n) < p). (4.1) 

If u is a local minimizer, under suitable regularity conditions for W, formally we obtain the 
following variation formula: 

f {W v (u(x))v(x) + V ( W(u(x)) ■ Vv(x)}dx = ( v t> G F C i?d(Q)). (4.2) 

This implies 

-div [V f W(«(x))] + W,,(u(x)) = in O. 

For fixed O and F C ff^O) as Problem HH we consider a family of minimization problem 
parametrized by tp G O(fi), where O(O) is defined by (|3.4jl . We define an affine space in 

Vfag) :=< P *(V(g)) = {veH 1 (n(v)); <p>\v) - g eV) (^eO(fi), g e H\n)). 

Problem 4.2. For given cp G O(Q) and g G if (fi), /md a local minimizer u(ip) of E(-,Q(tp)) 
in V(ip,g), i.e. u((p) G V(ip,g) and there exists p > such that 

E(u((p),n(<p)) <E(w,Q(tp)) ( ¥ we%fl) with Hw-uteOIUicnfe.)) </>)• (4-3) 
We define 

E*{ip) := E(u(<p),Q(ip)), (4.4) 
for a local minimizer it* (<£>). Using the formulas (|3.2p and (|3.3p . we define 

:= / W(<p(x), v(x), [A(ip)(x)]Vv(x)) K(ip)(x)dx (v G if 1 (0), peO(fi)), (4.5) 

then we have 

E(tp*v t Sl(tp)) = £{v,ip) (veH^Q)). 
The following theorem is a direct consequence of Theorem 12.11 

Theorem 4.3. Suppose that £ defined by (j4.5j) belongs to C°(ii 1 (0) x and i/iai £ (u, •) G 

C 1 (C(^)) M« G if 1 ^) and d v S G (^(if 1 ^) x 0(fi), (W 1 > DO (flo ) K Tl )) / )> where d v £ denotes 
the partial Frechet derivative of £ with respect to ip G Co C W 1 '°°(Qo,W ri ). Let Uq C V and 
Oq G be open subsets with ipo G Co, and we define U(ip,g) := <p*(Uo + g) C V(p,g). If 

n(-) G C°(Oo, H 1 ^)) and u((p) is a global minimizer of E(-,Q(tp)) in U((p,g) for all ip G Oq, 
then we have E* G C (Oq) and 

K(<p) = d^£{p- l u{<p),y). 

In particular, we have 

K(ip ) = d v £{u{jpQ)^o). 
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Proof. Let v(ip) := (p~ l u((p) — g G V and £(v, (p) := £(t> + g, cp). Then 14(99) is a local minimizer 
of f2(</?)) in V(<p,g), if and only if v{(p) is a local minimizer of £(-,(p) in V. We apply 
Theorem 12. II to £ and obtain I?* G C 1 (C'o) and 

£*M = d<p£(v(<p),(p) = d v S(v((p) +g,ip) = d v £ (<^~ 1 «(cp), <p). 

□ 

The minimizer u(ipo) is denoted by u hereafter. Under the suitable regularity conditions for 
W(£,7/,C), for ^ G W 1 ' 00 ^,^), we have 



c^ffa, tpo)\p] (4.6) 



— / W(x + tfi(x), u(x), [A(<po + tn)(x)]Vu(x)) K((po + tfj,)(x)dx 



n 



t=o 



= [ (V^(u) • n - (V c W(u)) T (X7fi T )Vu + lU(u)div^) dx. (4.7) 

Using the above formula, we also consider the inner variation, which is another type of 
variation for Problem 14.11 as follows. 



Theorem 4.4. Under the conditions of Theorem \4-3\ if u = u((po) is a local minimizer of 
ProblemU.lY then we have 



d^£{u,ip )[fi] =0 (p€ W 1 ' 00 ^), suppQi) C SI). 

Proof. We suppose that fi G W 1 ' oc (Sl) with supp(,u) C SI. We define (pit) = (po + tfj, for tel. 
From Proposition 13.21 and 13.41 if |^|Lip,f2 < 1> the corresponding pushforward operator (p(t)* is 
a linear topological isomorphism from H 1 (SI) onto itself. Moreover, from the formulas (|3.2p and 

a, 

lim\\ip(t)*u-u\\ H i (n) = 0, 

holds. Since 

E((p(t)*u,Sl) > E(u,Sl) = E((p(0%u,Sl) (|*| < |^|£j, >no ), 

we obtain 

^£(«,^ )[m]- 



0= ^E((p(t),u,Sl) 



= ■%:£(u,<p(t)) 
t=o az 



t=0 

□ 



5 Application to linear elliptic problems 

In this section, we consider a second order linear elliptic equation with a quadratic potential 
energy including the anti-plane displacement model of two dimensional linear elasticity. 

Example 5.1. We consider the following potential energy. 

w(H, 7?, C) = \ {CB(0C + KOv 2 ) - f(0v, (5.1) 

E(v,Sl) = J ji (V T vB{x)Vv + b(x)v 2 ) -f(x)vjdx, 
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where k E N U {0} and -B(£) is an n x n symmetric matrix (which is denoted by K™ym) and it 
satisfies 

B E C k (tt ,KyZ) and 3 (3 > s.t. ( T B(t)( > /3 |C| 2 ( V £ G n , V C G 

and 

bec k {n ,R), f e w k ' 2 (n ,R), (5.2) 

with the condition > for £ E f^o- We remark that 

V^,r?,C) = ^ \J2 V6 ii (e)CiCi + V6(e)r/ 2 j -V/(e)r/, 

w^.O = H0v-f(0, 

We suppose that Fp is a nonempty Lipschitz portion of d£l and that a bounded trace operator 
70 : H l (Q) L 2 (r D ) is defined and 

f:={«e 7o(«)=0}, 

is not empty. Then the minimization problem 14,31 corresponds to the following linear elliptic 
boundary value problem. 

-dw(B(x)Vu) + b(x)u = f(x) in Q(<p), 

u = g on Tfl, (BVm) • v = on<9il\r£>, 

where ^ is a unit normal vector on dVl. Prom the Poincare inequality and the Lax-Milgram 
theorem, there exists a unique global minimizer u{ip) G V(g,<p) for ip E 0(f2). 

Example 5.2. We consider the following potential energy which corresponds to the Poisson 
equation. 

W(Z,v,0 = l\(\ 2 -f(0ri, (5.3) 
E(v,n) = J(^\Vv\ 2 - f(x)v) dx, 
where k E N U {0} and / E W k ' 2 (tt , R). This is a special case of Example 15. 11 We remark that 
V*W(£, r/, C) = -V/(£)»7, r/, C) = -/(£)> V C W(£, t?, C) = C- 

Under the same boundary condition, the minimization problem 14.31 corresponds to the following 
boundary value problem of the Poisson equation. 

—Ait = f(x) in Sl(tp), 
du 

u = g on Fd, —— = on 9fi \ Fjj. 
ov 

In the case of n = 2, this represents the anti-plane displacement model of the isotropic linear 
elasticity. 
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Lemma 5.3. For the potential energy of Example \5.1i £ G C k (H l (£l) x 0(0)) and (|4.7j) holds. 
Proof. In the case of Example 15.11 £ becomes 

S(v,<p) = n±(A(<p)Vv) T B{i P {x)){A{<p)Vv) + ] i b(v(x))v 2 -f( V {x))v}^)dx 

= \^>x{v,y){x) + ^f 2 (v,<p)(x) - * 3 (v,<p)(x)J K(<p)(x)dx 

where 

<p) := (A{<p)Vv) T (B o cp) (A(ip)Vv) , V 2 (v, if ) := {b o (p)v 2 , * 3 («, := (/ o 

Under the assumptions, from Theorem 13.31 we obtain 

[(v,<p) i ^ A(v9)Vv] G C^-H^O) x 0(O),L 2 (O) n ), 

[^h^So^] g C fc (0(fi) > L 00 (n,R nxn )), 

i — ► w T Bw] G C rec (L 2 (fi) n x L°°(fl,R nXn ),L 1 (O)). 

From these regularities, it follows that \&i G C fc (ff 1 (r2) x 0(0), L-'-(O)). Similarly, from the 
following regularities 

[<p^boip}e fc (0(O),L°°(O)), [(b,v) ^ bv 2 } G C 00 ^ 00 ^) x ^(O^L^O)), 

[v^fopje c k (o(n),L 2 (n)), [(f,v) -> fv] g c°°(L 2 (n) x H\n),L\si)), 

#2 and ^ 3 also belong to C k {H l (£l) x 0(0), Since « G C°°(0(O), L°°(0)) from Theo- 

rem[331 we conclude that £ G C k {H l (9.) x 0(0)). □ 
We note that if O satisfies the cone property, the regularity conditions (I5.2p can be weakened 
by using the Sobolev imbedding theorem (see [TJ etc.): H l {0,) is continuously imbedded in 
L p (0), where p G [1, oo) if n = 2 and p G [1, 2n/ (n — 2)] if p > 3. The relaxed conditions are 

3 g>l, b G VF M (O ,lR), / G W fc ' 9 (O ,M), if n = 2, (5.4) 

in j.. 2n 

6 G W 'a (Oo, M), / G W (Oo,K), if n > 3. (5.5) 
More precisely, we have the following proposition. 

Proposition 5.4. IFe suppose that O satisfies the cone property. For the potential energy of 
Example\5J\with the condition ((53} or ((53]) in stead of ((Oft . £ G C fc (iJ 1 (0) x 0(0)) and KTb 
holds. 

Proof. We suppose n = 2 and the condition (|5.4p . For the q > 1 in (|5.4p we define := 
(1 — l/(/) -1 . From 

[^6 op] g fc (0(O), L 9 (0)), [«h»] g C 00 (-ff 1 (0),L 29 *(0)), 

[(6,u) ^ £w 2 ] G C°°(L 9 (0) x L^*(0),L x (n)), 
^ 2 G C fc (IT 1 (n) x O(0),L 1 (0)) follows. In the same way, from 

[</? i— > / o (^] g C fc (O(0),L 9 (0)), [v i-> u] g C^-H^O),!/"^)), 
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[CM -> fv] E C°°(L«(fi) x L'^L 1 ^)), 

^ 3 G ^(H 1 ^) x O(^),L 1 (0)) follows. 

Next we suppose n > 3 and the condition (|5.5|) . Put p := n/(n — 2), p* := (1 — = n/2, 

and (2p)* := (1 - l/(2p))- 1 = 2n/(n + 2). Then we have 

[v^v] g c^Cfr^nj.x^cn)), 

[y? i— > 6 o e C fe (0(ft),L p *(ft)), [(&,«) ^ 6w 2 ] G C°°(L p *(ft) x L 2p (ft), L 1 ^)), 

~ /op] G C fc (0(ft),L^*(ft)), [(/,*) ^ / V ] G C°°(L^)*(0) x L^L 1 ^))- 

Hence, from these regularities, we obtain *3 G C k (H 1 (fl) x 0(Q), L 1 (CI)). □ 
We state our results under the condition (15. 2|) hereafter. But we remark that the following 
results are valid even under the condition (|5.4p or (|5.5p with the cone property. We obtain the 
following theorem from Lemma 15.31 and Theorem 12.61 

Theorem 5.5. For the potential energy of Example \5.1\ E* G C k (0(Q)) holds, and if k > 1 we 
have 

E'^ )[n] = [ (V$W(u) ■ fx - (V f VF(u)) T (V/x T )Vn + W(u)divfj,) dx, 
Jn 

for // G W 1,co (CIq, M. n ), where u is the global minimizer of E(-,fl) in V(g). In particular, in the 
case of Example \5.Sl we have 

KifoM = J q {-(V/ • H)u - (V T u)(Vfi T )(Vu) + Q|Vn| 2 - /«) divpj dx, 

for fie W^°°(no,R n ). 

We consider the case of Example 15. H The global minimizer u belongs to H 1 (fl), but not to 
H 2 (fl) in general. The following boundary integral (|5.6|) is called J-integral etc.). 

Theorem 5.6. Under the assumptions in Example 1 5. 1\ with k = 1, if there exists a sequence of 
subdomains in which the Gauss- Green formula holds and 

oo 

fll C Cl-2 C • • • fl with \^jCli = Q, 
l=i 

and if the global minimizer u belongs to H 2 (Qi) for each I G N, then we have 

K(<PoM = Km / {W(u) n-v- (V C W(u) ■ v) (Vn • ^dH^ 1 , (5.6) 

for n G W 1,co (flo, R ra ), where v denotes the outward unit normal of dfli and 7Y n_1 denotes the 
n — 1 dimensional Hausdorff measure. 

Proof. Under the conditions, we have 

-div [V ( W(u)] + W v (u) = in L 2 (^), 

which is equivalent to 

-div(BVit) + bu = f in L 2 (ft,). 
14 



Let us denote the Hessian matrix of u by V 2 u. Since 



W(u{x))divn(x)dx = / W(u(x))n(x) ■ vdU n - v - / V x [W{u(x))] ■ fi(x)dx, 
V x [W(u{x))] = V^W{u(x)) + W v (u(x))Vu(x) + [V 2 u(x)]V c W(u(x)), 

Wi,(u(a;))Vw(x) • fj,(x)dx = / div [VfW(n)] Vu(x) • /j,(x)dx 

Jn, 

(V c W(u) ■ v) (Vu(x) • n(x)) dU n ~ x 
VjW(w)V (Vu(x) • dx 



n-1 



(y c W(u) ■ v) (V«(x) ■ n{x)) dUl 
V}W{u) {(V 2 u(x))fi{x) + (Vfi T )Vu} dx, 



V x \W(u(x))]-n(x)dx = / {V^W(u(x)) + [V 2 u(x)}V ( W(u(x))} • n(x)da 

Oi JQi 



(V c W{u) ■ v) (V«(x) • ^(x)) dHj 



n— 1 



V^(u) {(V 2 u(a;))/i(a;) + (V^ T )Vu} dx 

i 

(V c W(u) ■ v) (Vu(x) ■ n(x)) dH n ~ l 
+ [ {V^W(u)- fidx-VJW(u)(\7fi T )Vu}dx, 

we obtain, for ^ e W 1,DO (n 0) R"), 

EZ((po)\ji] = [ {V^W{u) ■ n- (V ( W(u)) T (X7fi T )Vu + W(u)divfi)dx 
Jn 

= [ (V f W(u) • A* - (V c Ty(u)) T (V / u T )Vu + Ty(u)div/i) dx 

+ [ (v e w(u) -/i- (v c iy(u)) T (v / u T )Vu)dx 

+ / W(u) • v dUlr 1 - [ V x [W{u)} ■ iidx 
= [ (V e W(u) • /i - (V c W(u)) T (V / u T )Vn + VF(u)div^) dx 
+ [ {W(u) ii ■ v - (X7 ( W(u) ■ v) (Vu ■ fJ.)}dH a 



in— 1 
'■x 



The first term tends to if I — ► oo. Hence, we have the formula (|5.6p . □ 
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In particular, for the case of Example 15.21 we have 

KM\jA = fonj^ |Q|V«| 2 -/«) ii-v-{Vu-u){V U -^dHT\ 

Theorem 5.7. Let Q be a bounded domain with C 2 ' -boundary. Under the condition of Exam- 
ple \5.2\ we assume that Tjj = dQ with g = 0. Then the following formula holds. 

K(<PoM = ~f \Vu\ 2 iM-vdHT\ 

1 JdQ 

Proof. From the regularity theorem of elliptic boundary value problems, we have u € H 2 (£2). 
So, we can choose fli = f2. We remark that 

Vu • (j, = (Vu • ■ v), |Vu • v\ = |Vu| on d£l, 

since u and its tangential derivatives vanish on the boundary. Hence we have 




This theorem can be generalized as follows. 
Theorem 5.8. Under the condition of Example \5.SX we assume that \i G W 1 ' oo (n ,M n ) satisfies 
supp(/i) n supp( g |p^-) = 0, suppO) fl T D n dQ \ F D = 0, 

3 : an open set ofW 1 s.t. O D supp(/u) and d£l PiO is C 2 - class. 
Then the following formula holds. 

K(<Po)\l4 = ~\ I |Vn| 2 /i • v dU n x - 1 - [ fun-v dU n ~\ 
1 Jv D Jan\r D 
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